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Abstract 

Certain complex-contour (a.k.a. quantum-toboggan) generalizations of Schrodinger's 
bound-state problem are reviewed and studied in detail. Our key message is that the 
practical numerical solution of these atypical eigenvalue problems may perceivably 
be facilitated via an appropriate complex change of variables which maps their multi- 
sheeted complex domain of definition to a suitable single-sheeted complex plane. 
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1 Introduction 



One-dimensional Schrodinger equation for bound states 



2m dx'^ 



Ipnix) + V{x) Ipnix) = En tpnix) , 



(1) 



belongs among the most friendly phenomenological models in quantum mechanics 
[1]. For virtually all of the reasonable phenomenological confining potentials V{x) 
the numerical treatment of this eigenvalue problem remains entirely routine. 

During certain recent numerical experiments [2] it became clear that many stan- 
dard (e.g., Runge-Kutta [3J) computational methods may still encounter new chal- 
lenges when one follows the advice by Bender and Turbiner [1], by Buslaev and 
Grecchi [5], by Bender et al [B] or by Znojil [7] and when one replaces the most 
common real line of coordinates x G M in ordinary differential Eq. ([T]) by some 
less trivial complex contour of a; G C{s) which may be conveniently parametrized, 
whenever necessary, by a suitable real pseudocoordinate s G M, 



Temporarily, the scepticism has been suppressed by Weideman [8] who showed that 
many standard numerical algorithms may be reconfirmed to lead to reliable results 
even for many specific analytic samples of complex interactions V{x) giving real 
spectra via Eq. 

Unfortunately, the scepticism reemerged when we proposed, in Ref. to study 
the so called quantum toboggans characterized by the relaxation of the most common 
tacit assumption that the above-mentioned integration contours C{s) must always 
lie just inside a single complex plane TZq equipped by suitable cuts. Subsequently, 
the reemergence of certain numerical difficulties accompanying the evaluation of the 
spectra of quantum toboggans has been reported by Bila [9] and by Wessels [10]. 
Their empirical detection of the presence of instabilities in their numerical results 
may be recognized as one of the key motivations of our present considerations. 
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Figure 1: The central segment of the typical PT— symmetric double-circle tobogganic 



2 Illustrative tobogganic Schrodinger equations 
2.1 Assumptions 

Whenever the complex integration contour C{s) used in Eq. (j2]) becomes topologically 
nontrivial (cf. Figures [H - H] for illustration), it may be interpreted as connecting 
several sheets of the Riemann surface j^ii^^'-i^^heeted) supporting the general solution 
^{general) (^r^^ of the Underlying complex ordinary differential equation. It is well known 
that these solutions ^(aenerai) ^^-^ non-unique (i.e., two-parametric - cf. [9]). From 
the point of view of physics this means that they may be restricted by some suitable 
(i.e., typically, asymptotic [H [5]) boundary conditions (cf. also Ref. 0)- In what 
follows we shall assume that 

(Al) these general solutions ipide^n^rai) (^^^^ Yiyq on unbounded contours called "tobog- 
ganic" , with the name coined and with the details explained in Ref. [7] ; 

(A2) our particular choice of the tobogganic contours 



curve of X G C^^^\s) with winding parameter k = 3 in Eq. ( ITOl) . This curve is 
obtained as the image of the straight line oi z e C^^\s) aX e = 0.250. 




:obogganic) 




will be specified by certain multiindex q so that (^(*o''099amc)^^^ ^ C^^^(s); 
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Figure 2: An alternative version of the double-circle curve of Figure [T] obtained at 
the "almost maximal" e = ^(c^***^"') - 0.0005 (note that ^(c"*^™^ ^ 0.34062502). 



(A3) for the sake of brevity our attention may be restricted to the tobogganic models 
where the multiindices g are nontrivial but still not too complicated. For this 
reason we shall study just the subclass of the tobogganic models 



2m dx^ 



M^) G l2(c(^)) (3) 



containing, typically, potentials 



+ 



or 



vg^{x) = V(^jco){x) = ix' + 



[x-iy (x + i)2 

G G 



(x-i)2 {x + iy 



F > 1 



G'>1 



(4) 



(5) 



with two strong singularities inducing branch points in the wave functions. 

In this manner we shall have to deal with the two branch points = ±1 in 

^{general) ^^-^ ^ In the language of mathematics the obvious topological structure of 
the corresponding multi-sheeted Riemann surface 7^(»T*"'*««'ieeteQ!) "punctured" 



at X 



(BP) 

(±) 



±1. In the vicinity of these two "spikes" we shall assume the generic. 



"logarithmic" structure oi n^rnuiusheeted) _ 
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Figure 3: The extreme version of the double-circle curve C^^^\s) at e ^ ^{critical) ^ 
2.2 Winding descriptors q 

The multiindex q will be called "winding descriptor" in what follows. It will be used 
here in the form introduced in Ref. [12] where each curve C^^\s) has been assumed 
moving from its "left asymptotics" (where s ^ — 1) to a point which lies below one 
of the branch points = ±1. During the further increase of s one simply selects 

one of the following four alternative possibilities: 

• one moves counterclockwise around the left branch point a^(f)''' (this move is 
represented by the first letter L in the "word" q), 

• one moves counterclockwise around the right branch point (this move is 
represented by letter _R), 

• one moves clockwise around the left branch point a^(f)'^ (this move is repre- 
sented by letter Q or symbol L^^ = Q), 

• one moves clockwise around the right branch point x^(^^p (this move is repre- 
sented by letter P or symbol R^^ = P). 

In this manner we may compose the moves and characterize each contour by a word 
g composed of the sequence of letters selected from the four-letter alphabet R, L, Q 
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Figure 4: The change of topology at e ^ ^(critical) ^j^gj-^ j^q^ (fTUj) starts giving the 
single-circle tobogganic curves C^^^\s) at k = 3. 

and P. Once we add the requirement of PT— symmetry (i.e., of a left-right symmetry 
of contours) we arrive at the sequence of eligible words g of even length 2N. 

At = we may assign the empty symbol g = ^ ot g = to the one-parametric 
family of the straight lines of Ref . ^ , 

= s-ie, e>0. (6) 

Thus, one encounters precisely four possible arrangements of the descriptor, viz., 

g e {LR,L-^R~\RL,R-^L-^} , N = l (7) 

in the first nontrivial case. In the more complicated cases where > 1 it makes 
sense to re-express the requirement of PT— symmetry in the form of the string- 
decomposition g = 0,[J where the superscript ^ marks an ad hoc transposition, 
i.e., the reverse reading accompanied by the L ^ R interchange of symbols. Thus, 
besides the illustrative Eq. ([7]) we may immediately complement the first nontrivial 
list 

n G {L,L-\R,R-^} , N = 1, 
by its A^ = 2 descendant 

^LL, LR, RL, RR, L ^R,R ^L,LR ^,RL '^,L ^R ^,R ^,R ^R ^j- 

(8) 
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Figure 5: The fully developed version of the single-circle tobogganic curve C*^^^^(s) 
obtained at k = 3 and e = 0.400. 

etc. The four "missing" words LL~^ , L~^L , RR~^ and R~^R had to be omitted as 
trivial here because they cancel each other when interpreted as windings |12] . 

3 Rectifications 
3.1 Formula 

The core of our present message lies in the idea that the non-tobogganic straight lines 
([6]) may be mapped on their specific (called "rectifiable" ) tobogganic descendants. 
For this purpose one may use the following closed- form recipe of Ref . [12] , 

M:{z e C(°)(s)) ^ (x G C^'\s)) (9) 

where one defines 

x = -i^(i-z^Y-l. (10) 

This formula guarantees the VT symmetry of the resulting contour as well as the 
stability of the position of our pair of the branch points. Another consequence of 
this choice is that the negative imaginary axis of 2; = —i\z\ is mapped upon itself. 

Some purely numerical features of the mapping ffTOj) may be also checked via 
the freely available software of Ref. [T3j. On this empirical basis we shall demand 
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that the exponent k will be chosen here as an odd positive integer, k = 2M + 1, 
M = 1, 2, . . .. In this case the asymptotics of the resulting nontrivial tobogganic 
contours (with M ^ 0) will still parallel the k = 1 real line C^^\s) in the leading- 
order approximation. 

3.2 The sequences of critical points 

The inspection of Figures [2] and [3] and their comparison with Figures H] and [5] reveals 
that one should expect the emergence of sudden changes of the winding descriptors g 
during a smooth variation of the shift e > of the initial straight line of z introduced 
via Eq. ([6]). Formally we may set g = g{e) and mark the set of the corresponding 
points of changes of g{e) by the sub- and superscript in 

The quantitative analysis of these critical points is not difficult since it gets per- 
ceivably simplified by the graphical insight gained via Figures |2] - H] and via their 
appropriately selected more complicated descendants. Trial and error constructions 
enable us to formulate (and, subsequently, to prove) the very useful hypothesis that 
the transition between different descriptors g{e) always proceeds via the same mech- 
anism. Its essence is characterized by the confiuence and "fiip" of the curve at any 
j = 1, 2, . . . , M in 5 = At this point certain two branches of the curve 

C^^\s) touch and reconnect in the manner sampled by the transition from Figure H] 
to Figure m 

The key characteristics of this fiip is that it takes place in the origin so that we 
can determine the point x^^^^^^^"'''^ = which carries the obvious geometric meaning 
mediated by the complex mapping fllOp . Thus, the vanishing = is to 

be perceived as an image of some doublet of z = 2;^'^''***™') or, due to the left-right 
symmetry of the picture, as an image of a symmetric pair of the pseudocoordinates 



^(critical) _j_ 



^{critical) 



At any k, = 2M + 1 the latter observations reduce Eq. ([6]) to elementary relation 

l = {l + [i{s-ie)ry (11) 
which may be analyzed in the equivalent form of the following 2M -|- 1 independent 



relations 

g2.iW(2M+l) = 1 + + ^ ^ ^ ^2 _ ^2 ^ 2 is £ . (12) 

These relations numbered by m = 0, ±1, . . . , M may further be simplified via the 
two known elementary trigonometric real and non-negative constants A and B such 
that 

In terms of these constants we separate Eq. f|T2|) into it real and imaginary parts 
yielding the pair of relations 

s^-e^-A = 0, 2se = B. (13) 

As long as £ > we may restrict our attention to the non-negative s and eliminate 
s = B/{2e). The remaining quadratic equation 

By{2e)^-e^-A = 

finally leads to the following unique solution of the problem, 



e = ^ \I-A + + B^ . (14) 
v2 

This formula perfectly confirms the validity and precision of our illustrative graphical 
constructions. 



4 Samples of countours of complex coordinates 

For the most elementary toboggans characterized by the single branching point the 
winding descriptor q becomes trivial because it is being formed by the words in a 
one- letter alphabet. This means that all the information about windings degenerates 
just to the length of the word q represented by an (arbitrary) integer [Hj. Obviously, 
these models would be too trivial from our present point of view. 

In an opposite direction one could also contemplate tobogganic models where a 
larger number of branch points would have to be taken into account. An interest- 
ing series of exactly solvable models of this form may be found, e.g., in Ref. [T5] . 
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Figure 6: The quadruple- circle tobogganic curve of x G C^^^^^\s). With winding 
parameter k = 5 in Eq. ffTOj) this sample is obtained at £ = g^^'^'^***'^"^) _ 0.0005, i.e., 
just slightly below the first critical value of g^^'^^***™') ^ 0.21574990. 

Naturally, the study of all of these far reaching generalizations would still proceed 
along the lines which are tested here on the first nontrivial family characterized by 
the presence of the mere two branch points in il){x). 

From the pedagogical point of view the merits of the two-branch-point scenario 
involve not only the simplicity of the formulae (cf., e.g., Eq. ( |T0|) in preceding sec- 
tion) but also the feasibility and transparency of the graphical presentation of the 
integration contours C*^^-* of the tobogganic Schrodinger equations. This assertion 
may easily be supported by a few explicit illustrative pictures. 



4.1 Rectifiable tobogganic contours with k, = ?> 

The change of variables ( ITOll generating the rectifiable tobogganic Schrodinger equa- 
tions must be implemented with due care because the knot-shaped curves C^^\s) 
may happen to run quite close to the points of singularities at certain values of s. 
This is well illustrated by Figure [T] or, even better, by Figure El At the same time 
all our Figures clearly show that one can control the proximity to the singularities 
by means of the choice of the shift e of the (conventionally chosen) straight line of 
the auxiliary variable z G C-^^ given by Eq. (jH]). 
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Once we fix the distance e of the complex line C^*^-* from the real line M we may 
still vary the odd integers /t. Vice versa, even at the smallest k = 3 the recipe 
enables us to generate certain mutually non-equivalent tobogganic contours C^^\s) 
in the e— dependent manner. This confirms the existence of discontinuities. Their 
emergence and form are best illustrated by the pair of Figures [3] and HI 
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(^'''''''''''''''^ 


——-^^^^^ Re X 

0,5 ^ 











Figure 7: The topologically different, triple-circle curve C'^^^^^\s) obtained at = 5 
and e = ej'^"*-'^') + 0.0005. 

We may conclude that in general one has to deal here with the very high sen- 
sitivity of the results to the precision of the numerical input or to the precision of 
computer arithmetics. This confirms the expectations expressed in our older paper 
[12] where we emphasized that the descriptor q is not necessarily easilly inferred from 
a nontrivial, detailed analysis of the mapping M.. 

4.2 Rectifiable tobogganic contours with > 5 

Once we select the next odd integer k = 5 in Eq. (ITOl) the study of the knot-shaped 
structure of the resulting integration contours C*^^-'(s) becomes even more involved 
because in the generic case sampled by Figure [6] the size of the internal loops proves 
unexpectedly small in comparison. As a consequence, their very existence may in 
principle escape our attention. Thus, one might even mistakenly perceive the curve 
of Figure |6] as an inessential deformation of the curves in Figures [1] or [2j 
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Figure 8: The other extreme triple-circled k = 5 curve C^^^^^\s) as emerging at 
£ = _ 0.005, i.e., close to the second boundary g^^'^"**^"') ^ 0.49223343. 

Naturally, not all of the features of our toboganic integration contours will change 
during transition from k = 3 to k = 5. In particular, the partial parallelism between 
Figures [2] and E] survives as the similar global-shape partial parallelism between 
Figures H] (where k = 3) and [7] (where k = 5). Moreover, a certain local-shape 
partial parallelism may be also found between Figure |2] (where the two upwards- 
oriented loops almost touch at k = 3) and Figure M (where the two downwards- 
oriented "inner" loops almost touch at /t = 5). The latter parallels seem to sample 
certain more general mechanism since Figure H] also finds its replica inside the upper 
part of Figure [H etc. Obviously, the next-step transition from k = 5 to k = 7 (etc.) 
may be also expected to proceed along similar lines. 

For the computer-assisted drawing of the graphical representation of the curves 
C^^^ the formulae of paragraph 13.21 should be recalled as the source of the most useful 
information about the critical parameters. The extended-precision values of the 
underlying coordinates of the points of instability are needed in such an application. 
Their M < 6 sample is listed here in Table 1. 

On this basis we may summarize that at a generic k, the variation (i.e., in all of 
our examples, the growth) of the shift e makes certain subspirals of contours C-^) 
larger and moving closer and closer to each other. In this context our Table 1 could. 
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Figure 9: The twice-circling tobogganic k = 5 curve C^^^^^^{s) as emerging slightly 
above the second critical shift-parameter, viz., at e = ^^"^^^^'■^ _|_ 0.005. 

in principle, serve as a certain systematic guide towards a less intuitive classifica- 
tion of our present graphical pictures characterizing transitions between different 
winding descriptors g and, hence, between the topologically non-equivalent rectifi- 
able tobogganic contours C^^\ During such phase-transition-like processes |1] the 
value of e crosses a critical point beyond which the asymptotics of the contours are 
changing. As a consequence, also the spectra of the underlying tobogganic quantum 
bound-state Hamiltonians will get, in general, changed [16] . 

5 Conclusions 

We confirmed the viability of an innovated, "tobogganic" version of PT— symmetric 
Quantum Mechanics of bound states in models where the general solutions of the 
underlying ordinary differential Schrodinger equation exhibit two branch-point sin- 
gularities located, conveniently, at x^^^^ = ±1. 

In particular we clarified that many topologically complicated complex integra- 
tions contours which spiral around the branch points x^^^^ in various ways may 
be rectified. This means that one can apply an elementary change of variables 
z{s) —J- x{s) and replace the complicated original tobogganic quantum bound-state 
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problem by an equivalent simplified differential equation defined along the straight 
line of complex pseudocoordinates z — s — ie. 

In detail a few illustrative rectifications have been described where we succeeded 
in an assignment of the different winding descriptors q to the tobogganic contours 
controlled solely by the variation of the "initial" complex shift e. An interesting sup- 
plementary result of our present considerations may be also seen in the constructive 
demonstration of feasibility of an explicit description of these transitions between 
topologically non-equivalent quantum toboggans characterized by non-equivalent 
winding descriptors q. Still, the full understanding of these structures remains to be 
an open problem recommended to a deeper analysis in the nearest future. 

In summary we have to emphasize that our present rectification-mediated recon- 
struction of the ordinary-differential-equation representation of quantum toboggans 
could be perceived as an important step towards their rigorous mathematical anal- 
ysis and, in particular, towards the extension of the existing rigorous proofs of the 
reality/observability of the energy spectra to these promising innovative phenomeno- 
logical models. 
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Table 1: Transition parameters for n — 2M + 1 with M — 1,2, ... ,6 



M 


m 


B 


(critical) 

[critical shift in C^^\s)] 


pseu 


docoordi 

(critical) 


nate 


angle 

(critical) 


1 


1 


0.8660 


0.34062501931660664017 


1.2712 


0.2618 


2 


1 


0.9510 


0.49223342986833679823 


0.96606 


0.4712 




2 


0.5878 


0.21574989943840034163 


1.3622 


0.1571 


3 


1 


0.7818 


0.49560936234793313854 


0.78876 


0.5610 




2 


0.9749 


0.41300244005317039597 


1.1803 


0.3366 




3 


0.4339 


0.15634410200136762402 


1.3876 


0.1122 


4 


1 


0.6428 


0.47438630343334929661 


0.67749 


0.6109 




2 


0.9848 


0.47917814904271720218 


1.0276 


0.4363 




3 


0.8660 


0.34062501931660664017 


1.2712 


0.2618 




4 


0.3420 


0.12231697600600608108 


1.3981 


0.08727 


5 


1 


0.5406 


0.44984366535166445772 


0.60092 


0.6426 




2 


0.9096 


0.49834558687374848153 


0.91265 


0.4998 




3 


0.9898 


0.42964189183273983152 


1.1519 


0.3570 




4 


0.7557 


0.28670826353957054964 


1.3180 


0.2142 




5 


0.2817 


0.10037407570525388131 


1.4034 


0.071400 


6 


1 


0.4647 


0.42666576745054519911 


0.54460 


0.6646 




2 


0.8230 


0.49875399287559237235 


0.82504 


0.5437 




3 


0.9927 


0.47264256935707423545 


1.0502 


0.4229 




4 


0.9350 


0.38168235795277279438 


1.2249 


0.3021 




5 


0.6631 


0.24649719795540125795 


1.3451 


0.1812 




6 


0.2393 


0.085076232785825555735 


1.4065 


0.06042 
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